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We formulated the A(1405) (abbreviated as A*) — > (Svr)'' invariant-mass spectra produced in 
p + p— )-p + A*+ A'"*" reactions, in which both the incident channel for a quasi-bound K~p state 
and its decay process to (Stt)" were taken into account realistically. We calculated M(E7r) spectral 
shapes for various theoretical models for A*. They are asymmetric and skewed, and were compared 
with recent experimental data of HADES, yielding M(A*) = 1405l"MeV/c^ and T = 62±10 MeV. 
The nearly isotropic proton distribution observed in DISTO and HADES is ascribed to a short 
collision length in the production of A*, which justifies the high sticking mechanism of A* and the 
participating proton into K^pp. 

PACS numbers: 21.45.-v, 13.75.-n, 21.30.Fe, 21.90.-ff 
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I. INTRODUCTION 

The A(1405) resonance discovered in 1961 [H (called 
herein A*) has strangeness S = -1, spin-parity = {\)^ 
and isospin /=0. It has been interpreted as a quasi- 
bound state of K~p embedded in the E -|- tt contin- 
uum since Dalitz-Tuan's original prediction 01 ■ In recent 
years, Akaishi et al. derived phenomenologically a com 
plex KN interaction (called here A Y interaction) 
based on the mass and width of A(1405), M 
MeV/c2 and r = 50 ± 2 MeV (so-called A(1405) 

ansatz). They applied this very attractive interaction 
to few-nucleon systems involving one and two K's^ and 
found nuclear bound states with unusually high nuclear 
density 0, On the other hand, a totally different 

framework with a double-pole structure of A(1405) has 
emerged on the basis of chiral SU(3) dynamics (called 
here Chiral), which claims that A(1405) consists of two 
poles around 1420 and 1390 MeV/c^, which are coupled 
mainly to KN and Ett channels, respectively (Tsl. Il4|. 
Then, the resulting weakly attractive KN interaction 
leads to much shallower K bound states [l^ [11] . 

Thus, it is vitally important to determine the location 
of the K~p resonance, "whether A(1405) is located at 
1405 MeV/c^ or above 1420 MeV/c^", from experimen- 
tal data without prejudice. For this purpose we have 
to treat the A(1405) structure with the model and 
the Chiral model on equal footing to be compared with 
experimental data. To solve this issue, observations of 
M(E7r) spectra associated with resonant formation of A* 
in the stopped- if ~ absorption in ^■"'He [l^, and also in 
d [l^ have been proposed. Whereas old bubble-chamber 
experiments of stopped K" in ^He [l^ indicated prefer- 
ence of A(1405) over A(1420) [l,[ll], a much more precise 
experiment with a deuteron tar get i s expected at J- PARC 
|20| . Alternatively, Jido et al. [2l[ proposed an in- flight 
K~ reaction on d, whereas Miyagawa and Haidenbauer 
[23 | questioned about the effectiveness of this method. 



In any case, old data on the in-flight K^ -i- d reaction 
by Braun et al. [l^l had a large statistical uncertainty 
in distinguishing A(1420) and A(1405), according to our 
statistical analysis. Future experiments at J-PARC of 
both stopped- if" [l^l and in-flight K~ [23] on d are ex- 
pected to give a convincing conclusion. 

Recent experiments on high-energy pp collisions have 
produced important data on the production of A(1405): 



p + p^p + K* + 



A* 
A* 



(I.l) 



The ANKE experiment at COSY with an incident ki- 



netic energy (Tp) of 2.83 GeV by Zychor et al. [25 1 



has yielded a (E^tt")'^ invariant-mass spectrum. It was 
analyzed by Ceng and Oset [1^ based on chiral SU(3) 
dynamics. They showed that the reaction in the A* 
production region is dominated by the |T2ipg2 process, 
and claimed that the spectrum develops a pronounced 
strength around 1420 MeV/c^, which differs from the 



1405 MeV/c^ peak in Hemingway's data [27[ analyzed by 
the |T'22p(;2 process [1,0] (see also Akaishi et al. ps}). 
This result might have been accepted as evidence for a 
double-pole structure of A* predicted by chiral SU(3) dy- 
namics [H, [13], if the statistics of the data were good 
enough. The ANKE data was also analyzed by Esmaili 
et al. [l8| , who, on the contrary, showed from a fair statis- 
tical comparison between the two models that the data 
was in more favor of the A Y model, but the statistical 
significance was not sufficient to conclusively distinguish 
between Chiral and AY. Thus, new data from HADES of 
GSI, which has just been published [1^, 113], is valuable 
for solving the present controversy. 

In the present paper we formulate the spectral shape 
of the (Ett)" mass to provide theoretical guides to ana- 
lyze experimental data of (Ett)" mass spectra from the 
above reaction. We take into account both the forma- 
tion and the decay processes of A(1405) in pp reactions 
realistically, following our KN — Ett coupled-channel for- 
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malism method Q. In this way, we derive the general 
form of the spectral function, which is not symmetric, 
but skewed with respect to the pole position. Then, we 
analyze (I]+-7r-+)° spectra from HADES at Tp = 3.50 
GeV 120]. 



II. FORMULATION 




A. Coupled-channel treatment of A* 

Our coupled-channel treatment of A(1405) is described 
in 0, [l^. The transition matrices of the two coupled 
channels, ^A^(l) and S7r(2), obeys the following equa- 
tion with a loop function Gi : 



Tij — Uij + ^ UiiGiTij, 



(II.2) 



where (i,j)=(l,2) and (2,1). The solution is given in a 
matrix form by 



(11.3) 



In our treatment, the loop function is considered to be 



(II.4) 




where kj is the relative momentum in channel j, and s^j's 
are non-dimensional strength parameters [5] . 

Among the matrix elements, Tn, T12, T21 and T22, the 
experimentally observable quantities below the K + N 
threshold are -(l/7r)Im Tn, | T21 P 92 and | T22 P 92, 
where the second term is a Stt invariant-mass spectrum 
from the conversion process, KN Stt (we call this 
"T21 invariant mass"). The T21 invariant mass coincides 
with the KN missing-mass spectrum in the mass region 
below the K + N threshold, as denoted by relation [18j . 
that 



ImTii = iTairimGz. 



(11.5) 



The third term is a Ett invariant-mass spectrum from the 
scattering process, Ett — > Ett (we call this "r22 invariant 
mass" ) . 



B. A* (Stt)" spectrum shape 

The diagram for the reaction eq. (II. 1[) is shown in 
Fig.[TJ The decay processes via T21 and T22 are also given 
in this figure. The kinematical variables in the cm. of 
the p + p collision for both the formation and the decay 
processes are given in Fig. [H 

In the present reaction we use |T2ipg2 because the in- 
cident channel to bring A(1405) is K~ +p together with 
(see Fig. [J (b)). This was also concluded by Ceng 
and Oset l26|, who studied the reaction mechanism in 




FIG. 1: (Color online) Feynman diagrams for the p + p — >■ 
p + -I- A* — >■ p + + (Stt)'^ reaction: (a) for the process 
via T22 and (b) for the process via T2i. 



detail. The |T22p(72 spectrum would be applicable when 
E and tt mesons are available in the incident channel, 
as shown in Fig. [T] (a). The IT22P92 spectrum is char- 
acterized by a large tail [3| in the higher-mass region 
up to the kinematical limit, which can in principle be 
recognizable by an observed spectrum. Experimentally, 
however, a bump in the upper-tail region may be masked 
by an ambiguous shape of the continuous background, 
and may thus be difficult to extract. We may allow a 
small admixture of |T22p'72 in our likelihood analysis of 
the experimental data. 

The |r2ipg2 and |r22p92 curves of the Chiral model, 
as given by Hyodo and Weise [lB| as well as those of the 
model, are shown in Fig. 1 (upper) of Ref. [l^. They 
will be compared with the new HADES data at the end 
of the present paper. 



C. Spectral function in pp reaction: S{x) 

Now, we consider the spectrum function of the invari- 
ant mass, S{x), in the case oi pp reactions. We compose 
it by the incident channel, Wform(x), and the decay chan- 
nel, G{x), as follows: 



with 



S{x) = VKfo™(x) X G(x), 



X = Af(E7r). 



(H.I 



(11.7) 



G{x) is expressed in terms of the T matrices, T22 and 
T21, as shown in Fig.[T](a), (b). Each function calculated 
for an assumed M of the A* pole is shown in Fig. [31 
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FIG. 2: (Color online) Kinematical variables in the center of 
mass of the p+p collision for (a) the formation process, VKform, 
and (b) the decay channel, G{x). 



D. Formation process function: H^form 

The A* formation probability in a pp collision using 
the T-matrix is given as follows: 



W^form(a;) 



2 Tn 



Eo 



(27r)3(fe)6 ko 



X hk2EiE2 [1 + 1^(1 + ^ cos(V+)] 

where i?o and fcg are the initial energy and momentum 
in the c. m. frame, as given by 



fci 



fco 



(11.9) 



Also, Opx+ = (^p ~ (^K+) is the angle between K'^ and p, 
6 is the range of the p — p reaction, and the momentum 
transfer, Q, is 



Q = [fco^ + k2^ -2 ko k2 cos 9p]2. 



(11.10) 



As can be seen from the factor, 1/(1 + b^Q^)^, a shorter 
range of b can effectively moderate the strong suppression 
due to a large momentum transfer, Q, in a high-energy 
pp collision. 

Figure [3] (b) shows the behavior of Wform(2;) for Tp = 
2.50, 2.83 and 3.50 GeV, which are normalized at a; = 



2.5GeV 
2.8GeV 
3.5GeV 
/W=1405 MeV/c 




—I 1 1 I 1 1 1 1 1 1— I 1 1 1 1 1 1 1 1— 

1300 1350 1400 1450 1500 
M(I.n) [MeVIc] 

FIG. 3: (Color online) Normalized spectral functions S{x) (a) 
composed of the formation-process function W^form (b) and the 
decay-process function G{x) (c) for Tp — 2.50, 2.85 and 3.50 
GeV. TUB = 770 MeV/c^ and {ep,ej,K+) = (90,180). The M 
value of A* is assumed to be 1405 MeV/c^, as indicated by 
the vertical dashed line. 



1400 MeV/c^. They have respective kinematical upper 
limits, and make the mass distribution damp toward the 
kinematical limit. As a result, the observed spectrum 
shape, S{x), changes, as demonstrated in Fig. [3] (a), 
whereas G{x) is independent of Tp. 



E. Decay process function: G(x) 

The decay rate of A(1405) to (Stt)" is calculated by 
taking into account the emitted S and tt particles real- 
istically, following the generalized optical potential for- 
malism in Feshbach theory [s^l, given by Akaishi et 
al. 0, [2^. The decay function, G{x), is not simply a 
Lorentzian, but is skewed because the kinematic freedom 
of the decay particles is limited, particularly, when the 
incident proton energy, Tp, decreases and approaches the 



4 





1 i i 1 1 i i i 1 1 

7=3.5 

p ^ 


' ' ' ' 1 ' ' 
3eV {xO.5) 






M=1407 MeV/c' 




/ 3.0 








"c 






// 














55" 


// / 






/^^-"^^^^-^ 2.7^ 






^ — 






^^^^^^^^^^ ^ 2^ 






1 1 ^ 1 1 1 1 1 1 1 


"1--^^^ ■ — 1 — ■ — r-l 



1300 



1350 1400 

MM(pK*) [MeV/c^] 



1450 



FIG. 4: (Color online) Incident energy dependence of the ab- 
solute values of the spectral function at ms ~ 770 MeV/c^ 
and (^p, V+) = (90,180). 
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production threshold. Its general form is given as 



X I (fc(a;)k|fco(a;)) I , 



(11.11) 



where the relative momenta in the entrance and exit 
channels of Fig. [5] (b) are calculated by: 



and 



with 



ko{x) 



~k{x) 



c^X{x, rriK, Mp 
2hx 



2hx 



(11.12) 



(11.13) 



X{x, mi, 7712) — {x + mi + m2){x + mi — ni2) 



X [x ~ mi + m2){x — mi — m2). 



(11.14) 



It should be noticed that A(a;, mK,Mp) becomes nega- 
tive at around x = 1400 MeV/c^, where we must choose a 
positive Im k on the physical Riemann sheet. This case 
corresponds to direct excitation of the A* quasi-bound 
state from the p + p channel. 

In the case of AY, the T— matrix is 



(k\t2i\ko) = g{k) T21 g{ko) 
for the T21 channel, and 

9ik) 



A2 +fc2 



(11.15) 



(11.16) 



with A = mBc/h, ms being the mass of an exchanged 
boson, and k is the relative momentum of S and tt. 



FIG. 5: (Color online) Experimental spectra of AM{pK^) 
in the pp -> pAK^ reaction at Tp = 2.50 and 2.85 GeV in 
DISTO experiments. Taken from [33] • 



The shape of G{x), as given by eq. ljll.lip . includes the 
momenta, feg and k, which are functions of Tp. How- 
ever, the function G{x) is shown to depend only on the 
invariant-mass x, namely, G{x) is a unique function of x 
and does not depend on Tp. It is bounded by the lower 
end {Ml = m^) = 1328 MeV/c^) and the upper 

end {Mu = Mp + mK-) = 1432 MeV/c^). 

It is to be noted that the position of the peak in G{x) 
is significantly lower than the position of the pole {M — 
1405 MeV/c^) in T21, as assumed here and indicated by 
the vertical dashed line. Furthermore, the position of the 
peak (or centroid) of S [x) is lowered due to the formation 
channel function Wform(a^)- 



III. NUMERICAL RESULTS 

In this section we present results from numerical calcu- 
lations, and discuss their physical implications. The im- 
portance of the present work is to consider both Wforai{x) 
and G{x) functions. In most illustrative samples, we ap- 
plied the A Y model with the PDG parameters of [3] , 
=1407 MeV/c^ and F = 50 MeV. To compare the Chi- 
ral model with the yiy model on equal footing, we also 
applied the same procedure as above to Hyodo-Weise's 
T-matrices to obtain realistic spectrum shapes iS'(a;). 



A. Dependence on the incident energy, Tp 



Regarding eqs. pi.61 III.8[ IIL11[) again, it is clear that 
the spectral function depends on the incident proton en- 
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ergy due to the Wform(a;) function and G{x). Figure |4] 
shows absolute values of spectral functions S{x) for var- 
ious incident energies (Tp) at tob = 770 MeV/c^ and 
{0p,OpK+) = (90,180). The shape of S{x) is nearly the 
same, but toward the reaction threshold (^Tp^'^°^^ = 2.42 
GeV) not only the absolute value diminishes, but also the 
spectral shape changes drastically, as shown in Fig. [Sja) 
for the normalized spectral functions at Tp =3.50, 2.83, 
2.50 GeV. The most extreme case is seen at Tp — 2.50 
GeV, where the main part of x > 1400 MeV/c^ is miss- 
ing due to the kinematical constraint, and a very skewed 
component below 1400 MeV/c^ appears. 



B. Behavior near the production threshold of Tp 

The above prediction is indeed in good agreement with 
the observed spectra of DISTO at Tp = 2.50 and 2.85 
GeV [33|, as shown in Fig. [5] Even in such a very skewed 
spectrum, one can extract the decay function, G{x), from 
an observed spectral function by taking the ratio 



DEV[G{x)] 



_ Six) 



obs 



(III.17) 



using a calculated Wform function. This is a kind of 
the deviation spectrum method introduced in stopped- iiT^ 
spectroscopy |18| . 



C. Angular distribution and correlation 

The cross section of this reaction has substantial an- 
gular dependence (Fig. [T]), but the bound-state peak is 
distinct at any angle, and we can choose (^p,^p_R-+) = 
(90,180), because the cross section is modest and the 
peak-to-background ratio remains large. The normal- 
ized cross sections (spectral shapes) at various angles are 
found to be nearly the same. Since the two incident pro- 
tons are indistinguishable, the A(1405) formation process 
is angular symmetric, as shown in Fig. [T] We can write 



JpK 



-) 



(III.18) 



for Op^Q^ 90° and_0^+ = - 180°. 

According to eq. pi.8p and eq. (jII.10|) . Wfomi, and thus 
the spectral function, S{x), are related to the outgoing 
proton angle. Op, and the angle between the outgoing 
proton and , 0pK+i as shown in Fig. [51 Although 
these curves look different, the spectrum shape does not 
depend on the angle. We choose and use 9p — 90°, 
9pK+ = 180° in all of the following calculations. 



D. Dependence on the exchanged boson mass 



c 
CO 
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FIG. 6: (Color online) The spectral functions for various an- 
gles, {Op, epji+), for Tp = 3.50GeV and ms=770 MeV/c^ 




15 30 45 60 75 90 105 120 135 150 165 180 



FIG. 7: (Color online) Normalized angular distributions of the 
outgoing proton for different exchanged boson masses, ms = 
2000, 770 and 140 MeV/c^ at Tp = 3.50GeV. 



Tp = 3.50GeV. The nearly isotropic angular distribution 
with a large boson mass explains the experimental data 
of HADES at Tp = 3.50 GeV [H, H^, which shows that 
the proton angular distributions together with A(1405) 
and A(1520) are nearly isotropic. A similar behavior is 
observed in the DISTO data at Tp = 2.85 GeV (see Fig.[S] 
[H, US)- Such a short collision length as revealed in the 
production of A(1405) in the pp reaction is one of the key 
mechanisms (A* doorway) responsible for forming K~pp 
from high sticking of A* and p [1]. On the other hand, 
it is well known that the proton emitted in the ordinary 
— > p + A + K'^ reaction has sharp forward and back- 
ward distributions, indicating that the mediating boson 
is niB — niTr 33 - 1341 . 



IV. 



In this section we analyze the recent HADES data for 



Figure [7] shows the normalized angular distributions charged final states of E 7r+ and E"'"7r in a pp collision 



of the outgoing proton, 6p, for various masses of the ex- 
changed boson, niB = 2000, 770 and 140 MeV/c^, at 



at Tp = 3.50 GeV. The data we use are the missing-mass 
spectra, MM{pK~^), deduced by the HADES group, as 
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FIG. 8: (Color online) Comparison of HADES data (E+tt" + E"7r+, closed squares) at Tp = 3.50 GeV l30[ with best-fit 
theoretical spectral functions S{x). a) Best-fit HKAY curves (with = 9.5, M = 14051" MeV/c^ and T = 62 ± 10 MeV). b) 
model with the PDG parameters (with = 14, M = 1405.lj:i;o MeV/c^ and T = 50 MeV ;8]). The Chiral model using 
HW's T21 (with = 111, c)) and T22 (with x^ = 40, d)). 
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FIG. 9: (Color online) Confidence level contours from fit- 
ting of the HADES data of E+vr" + E-7r+ at Tp = 3.50 GeV. 
The PDG values are also shown. 



given in Fig. 1 of [30| , which are corrected for acceptance 
and efficiency of the detector system. They are expressed 



as 



Y{x) = Fa- {x) + Fe- {x) + Yxiz2q{x) + Yno,,Kos{x), 

(IV. 19) 

with Fa* for A*, Fs- for S(1385), rAi520 for A(1520), 
and InohRcs for non-resonant continuum. They decom- 
posed the experimental data, Y{x)^ by the above four 
components, which were obtained by model simulations, 
among which the I](1385) and the A(1520) components 
were determined by using the experimental data. The 
shape of the non-resonant Stt continuum was simulated. 
In their fitting they cautiously excluded the area around 
1400 MeV/c^ for MM{pK+) in order not to bias the fi- 
nally extracted shape of the A* resonance. Then, they 
found that a simulation of the A* region by using a rela- 
tivistic s-wave Breit-Wigner distribution with a width of 
50 MeV/c^ and a pole mass of 1385 MeV/c^ can repro- 
duce the experimental data very well, but using instead 
the nominal mass of 1405 MeV/c^ fails. 

This conclusion depends on their assumption of the 
symmetric Breit-Wigner shape, which is not valid in 
the case of a broad resonance with adjacent endpoints, 
M(S -I- tt) and M{p + K~), as we have seen. Thus, in 
turn, we decided to set up an excess component, Y\t(x), 
by subtracting the given three components from the ex- 
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perimental spectrum Y{x) as 

Ya* {x) = Y{x) - Fs* (x) - yAi52o(a;) - >NonRos(a;), 

(IV.20) 

where the statistical errors of Y{x) are inherited to 
YA^ix). 

The HADES spectra, as given in Fig. 1 of [3^, indi- 
cate that the spectra of both charged channels are similar 
to each other, yielding nearly the same M values. This 
fact indicates that the Ett resonance is formed by nearly 
pure charged states, E+tt^ and I]^7r+, without isospin 
mixing. It also justifies the use of T21 for the analy- 
sis of Af(I]7r) in the case of pp reactions. On the other 
hand, the statistical fluctuation of each charged-channel 
spectrum is rather large. Thus, for the final analysis we 
use the sum data of HADES (S+7r~ -I- E~7r+), which is 
presented in Fig. 1(c) of 30:]. Keeping the last three 
components of eq. (jIV.20p fixed, we fit the experimental 
data of Ya'{x) with ri =21 data points in the range of 
1300 to 1550 MeV/c^ (closed points with error bars in 
Fig. HI) by assumed theoretical functions S{x). 

Generally, the experimental histogram, Ni, i— 1, n 
with respective statistical errors, cr^, is fitted to a theoret- 
ical curve, S{x; M, T) with x = MM{pK^) involving the 
mass M and width F as free parameters by minimizing 
the value: 

(=1 \ / 

Figure [5] shows the results of the fitting, where the 
HADES data (E+tt" +^-n+) at Tp = 3.50 GeV [s^ 
are compared with best-fit theoretical spectral func- 
tions, S{x). The present AY treatment (hereafter, called 
HKAY), adopting the PDG values (M = 1405. 1^^:^ 
MeV/c^ and G = 50 MeV [&]), gives a remarkable fitting 
with — 11, which is comparable with the statistically 
expected value, < >cxp~ 19. On the other hand, the 
Chiral model gives much larger x^ values of ~ 111, when 
T21 is chosen, and of 39, when T22 is chosen. Another 
Chiral model spectrum by Geng and Oset [1^ is almost 
identical to HW's T2i. Thus, the chiral models indicate 
a substantial deviation from the experimental data. 

Furthermore, we can find best-fit values of (M, F) from 
drawing confidence contour curves by varying the param- 
eters, (M, F) in a plane. The results are shown in Fig. |9l 
From this contour mapping we obtain the following best- 
fit values with 68% CL (Ict) errors: 

M = 1405i^^ MeV/c^ (IV.22) 
F = 62 ± 10 MeV. (IV.23) 



The best-fit curves are shown together with the experi- 
mental points in Fig. [5] The M value thus obtained from 
the present analysis of the new HADES data confirms 
the traditional value 0, [1] . 

V. CONCLUDING REMARKS 

We have presented results of our calculation for the 
spectral shape of MM{pK^) in the pp — ?> pA*K^ reac- 
tion based on the KN-YiTt coupled-channel treatment. 
We took into account both the entrance process and 
the decay process. The formation probability, VFtonm 
of A* in a pp collision and the decay rate, G{x) to 
(Ett)'^, were formulated. The spectral function is given 
by S{x) — Wfom X G{x). It was found to be asymmetric 
and skewed due to the kinematic limitation imposed by 
the entrance channel. The peak of S{x) is not located at 
the pole position. 

With this tool in hand we analyzed the recent HADES 
data. Although the observed spectra of M AI {pK^) ap- 
pear to show the peak position at around 1385 MeV/c^, 
the x^ fitting by our theoretical spectral functions pro- 
vided M = 1405^9^ MeV/c^. This value is in good agree- 
ment with the values obtained from a recent analysis Il7[ 
of an old experimental data of stopped-X" in ^He [ig |. 
taken up as the updated PDG value (M = lAOb.lt\:l) 
dl. On the other hand, the Chiral model with M ~ 1420 
MeV/c^ cannot reproduce the experimental data. 

The proton angular distribution in the A* production 
was also calculated. The isotropic distribution observed 
in HADES (s^ and DISTO [H, H were explained by 
a short-range collision with an intermediate boson mass 
heavier than the p meson mass. This is consistent with 
the calculated large cross section for the production of 
K~pp in pp collisions Q, which has recently been ob- 
served in DISTO experiments f32l|. 
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